Density matrix interpretation of solutions of Lie-Nambu equations 
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The spectrum of a density matrix p(i) is conserved by 
a Lie-Nambu dynamics if p(t) is a self-adjoint and Hilbert- 
Schmidt solution of a nonlinear triple-bracket equation. This 
generalizes to arbitrary separable (positive- and indefinite- 
metric) Hilbert spaces the previous result which was valid for 
finite-dimensional Hilbert spaces. 



I. STATE VECTORS VS. DENSITY MATRICES 
IN NONLINEAR QUANTUM MECHANICS 

There exist prejudices concerning nonlinear generaliza- 
tions of quantum mechanics. One of them is a belief that 
any generalization must lead to unphysical effects such 
as a faster-than-light transfer of information. Although 
the proofs of these unphysical phenomena are explicit 
and mathematically correct [Q-Q they are based on some 
physical assumptions which are unjustified. One of these 
physically wrong elements is a naive use of the projection 
postulate. Today we understand that if the dynamics is 
nonlinear we are not allowed to simply project a solution 
on some direction in a Hilbert space. One of the reasons 
is that a projection of a solution is in general not a solu- 
tion of a nonlinear Schrodinger equation. A more subtle 
argument is provided by the notion of nonlinear gauge 
transformations introduced by Doebner and Goldin pJq] 
and developed by the Clausthal school [Q . It is clear that 
there exists a class of nonlinear Schrodinger equations 
which arc obtained by a nonlinear gauge transformation 
from an ordinary linear Schrodinger equation. They not 
only give the same probability density in position space 
but also may look "truely" nonlinear (there exist nonlin- 
ear gauge transformations that simply add a nonlinear 
term but do not alter the kinetic and potential parts in 
a Hamiltonian) . Obviously if the nonexistence theorems 
were true such "nonlinear" equations would have to lead 
to unphysical effects. But the point is that they do not 
lead to any new effects since, by definition, they are phys- 
ically fully equivalent to to the linear theory. Therefore 
the nonexistence theorems must contain some elements 
which are physically wrong. From the perspective of the 
nonlinear gauge transformations it is clear that one of 
them is a wrong use of the projection postulate. Let us 
note, however, that the explicit example discussed in H 
is not based on this postulate (and thus is not equivalent 
to the examples given in [0J|]; a simple argument shows 
also that the "telegraph" discussed in II works in the op- 
posite direction than those from An element which 



is physically wrong here is the wrong way of describing 
composite systems. This was clarified by Polchinski || 
and Jordan j9|. The latter work was based on a density 
matrix reformulation of Weinberg's nonlinear quantum 
mechanics ]To| , |i~t| ] . 

Density matrices play in nonlinear quantum mechanics 
a role which is somewhat ambiguous. On the one hand, 
one of the earliest attempts of formulating a general non- 
linear framework for quantum mechanics was Mielnik's 
convex formalism" 121. Its main idea was to keep a 
figure of states convex and derive a probability interpre- 
tation in terms of its global geometric properties. From 
this perspective the density matrices might be even more 
fundamental than state vectors. On the other hand, how- 
ever, all works that start from pure states and nonlinear 
Schrodinger equations lead to difficulties when it comes 
to "mixtures" . The difficulties arc so deep that some au- 
thors tend to reject the very notion of a density matrix 
in a nonlinear context jL|, although different proposals 
of combining mixtures with nonlincarity of pure states 
exist in the literature (cf. |l^Jl5[| ). 

A nonlinear extension of quantum mechanics based 
on a triple bracket Nambu-type generalization of the 
Liouville-von Neumann equation (cf. |l6[|) proposed by 
one of us Jl7]-[l9| starts from a completely different per- 
spective. The idea is to find a general scheme which on 
one hand includes the linear and Weinberg-Bona-Jordan 
cases and on the other leaves some room for nonlinear 
generalizations that do not use nonlinear Hamiltonians. 
Such a starting point is motivated by ambiguities in prob- 
ability interpretation caused by the notion of eigenvalue 
of a nonlinear operator (20). The Lie-Nambu scheme 
proved very powerful and elegant and has, in our oppin- 
ion, several advantages over the standard paradigm of 
nonlinear Schrodinger equations. The density matrices 
play in this formalism a fundamental role. Still the basic 
question of an interpretation of solutions p(t) as density 
matrices was not fully clarified in the earlier work. The 
Theorem 5 discussed in Jl8[ worked essentially in finite di- 
mensional cases whereas the generic infinite dimensional 
Hilbert space problem was left open. In this paper we 
give an alternative proof which generalizes this theorem 
to infinite dimensional Hilbert spaces. We consider also 
the indefinite-metric case which is of some interest for a 
relativistic theory. 
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II. TRIPLE-BRACKET EQUATIONS AND 
SPECTRAL PROPERTIES OF THEIR 
SOLUTIONS 

Consider a one-parameter family p = p(t), t G R, of 
Hilbert-Schmidt self-adjoint operators acting in a separa- 
ble Hilbert space and satisfying the Lie-Nambu equation 



ipa = {p a ,H, S}. 



(1) 



Here p a := pAA'(d,(i') are components of p in some 
basis, A and A' are discrete (say, spinor) indices and 
a, a' the continuous ones (coresponding to, say, posi- 
tion or momentum). The dot represents a derivative 
with respect to the parameter t. In nonrelativistic case 
this is just an ordinary time. In the relativistic case 
the meaning of t depends on a formalism (t is time 
in some reference frame in [H, and a "proper time" 
in the off-shell formulation given in H = H{p) 

is any (functionally) differentiable functional of p and 
S = S{p) = S(C\{p), . . . , C k(p), ■ ■ ■) is differentiable in 
C n (p) (see Appendix V_B). We assume the following 



summation convention for the composite indices a: A 
contraction of two composite indices means simultane- 
ous summation over the discrete indices and integration 
(with respect to an appropriate measure) of the continu- 
ous ones. The triple bracket itself is defined as 



{F, G, H} = n c 



be 



5F 5G SH 

Sp a Spb Sp c 



(2) 



where O a j, c are structure constants of an infinite- 
dimensional Lie algebr a wh ich also depends on the model 
| p8| , pj[ (see Appendix VA). The indices in O a f, c can be 
raised and lowered by a metric discussed in detail in 
p8| , ^9| . The metric is well defined for both finite- and 
infinite-dimensional Lie algebras and is not equivalent to 
the Cartan-Killing one (the latter does not exist in the 
infinite-dimensional case). 

Before we proceed with the main theorem we shall first 
prove a few useful technical results. 

Lemma 1: Let {pk}kLi be a sequence of nonnegative 
numbers such that p\ > P2 > P3 > • ■ ■ and the series 
J2T=iPk i s convergent. Then 



lim 




Pi- 



Proof. Using the three sequences theorem one imme- 
diately generalizes the standard proof known for finite 
sequences. □ 

Lemma 2: Let {pk\kLn telfcLi be two sequences ful- 
filing the assumptions of the above lemma. Suppose that 
for every m £ N 



Then pk = qu for every k — 1,2,.... 

Proof. The equality p\ — q\ follows directly from Lemma 

1. The proof is completed by induction. □ 

Lemma 3: Let {qk}kLi, <7fc € R, be an arbitrary se- 
quence and {pk\uLi satisfy the assumptions of Lemma 1. 
Suppose that for every m 6 N 



k=l 



k=l 



Then {pk}"kLi = {?fe}feLi U P to permutation. 
Proof. Define two sequences {pk}kLn felfcLu where 
Pk = Pki 9k = ql- {<7fc}^Li is absolutely convergent so we 
can rearrange its terms in such a way that the rearranged 
sequence satisfies assumptions of Lemma 1. Assume this 
done. Lemma 2 implies that \qk\ — Pk, for any k. Let 
Q + be the subset of positive elements of {qk}k%i- Q+ 1S 
non-empty since otherwise all elements of {qk}^! wou ld 
be < which contradicts the assumption that J^ fc—1 qk — 
YlkLiPk > 0- (We exclude the trivial case where all pk — 
0.) Therefore for any qk G Q+ there exists pk such that 
Pk = qk- Subtracting these elements from both sides of 



fe=i fe=i 



we get for m = 1 



E 



Pk 



remaining pfc 



E 

q k <0 



qk 



which can hold if and only if all such = 0. Therefore 
all qk > and pk = qk for any fc.D 

We are interested in solutions of (Q) where p(t) are 
Hilbert-Schmidt self-adjoint operators acting in a sepa- 
rable Hilbert space. The following theorem states that 
spectrum of p(t) is conserved by the Lie-Nambu dynam- 
ics. 

Theorem 4 '■ Let p(t) be a Hilbert-Schmidt self-adjoint 
solution of (|l]) whose spectrum is sp p(t) = {^k(t)} k % 1 . 
If {Afc(0)}^L 1 satisfies assumptions of Lemma 1 then 
Afc(t) = Xk(0) for any t. 

Proof: For any t the solution can be written as p a (t) 

E 



k^iXk(t)(j) A (t,a)(j)\,(t,a'). According to Theorem 4 



in (18| the functional C„(p), n £ N, (see Appendix 
is time independent. Therefore 



E(^Tw i 

k=l 



fe=l 



(3) 



for any t and n (r / kk includes the indefinite metric case, 
see Appendix |VB| ) . Lemma 3 implies that the sequences 
{ Afc (i)}fe^i and {A/ c (0)}^° =1 are identical up to permuta- 
tion if the metric rj kl is positive definite. In the indefinite 
metric case we define A^(i) = Afc(i) 2 . ([|) implies 



Lemma 3 again implies that the sequences {|Afc(t)|}^_ 1 
and {Afe(0)}^ 1 are identical up to permutation. Finally 
continuity in t means that Afc(t) = Afe(O) in both cases. □ 



III. CONVEXITY PRINCIPLE: AN EXAMPLE 

A nonlinearly evolving density matrix cannot satisfy 
an ordinary convexity principle. Still, being a Hilbert- 
Schmidt operator it can be spectrally decomposed and 
the spectral projectors can be regarded as its "pure state 
components" . This is justified by the fact that the spec- 
trum of p(t) is time- independent. To sec how this works 
consider a simple example. The example simultaneously 
illustrates the peculiarity of the triple-bracket formalism: 
An existence of nonlinearities that become invisible on 
pure states. Let H (p) = Tr (hp) where h is a 2 x 2 Her- 
mitian matrix, p = /3q1 + p ■ cr, and 



of two not necessarily mutually orthogonal density matri- 
ces. Let pi(0) = p\q1 + pi • <t and P2(0) — P20I + P2 ■ f- 
The solution we look for is given by (Q) but with p = 
Pi Pi + P2p2- The Hilbert-Schmidt vectors <fr\ depend 
now not only on the eigenvalues of pi(0), p 2 (fy, but also 
on pi and p 2 . This implies also that the dynamics of the 
density matrix can be written here as 



p(t)= Pl Up 1 (0)tf + p 2 Up 2 (0)tf 



(8) 



where U — U(t, pi(0), p 2 (0),pi,p 2 ) is unitary but 
parametrized by the initial condition. 

Typically a nonlinear evolution is a result of a mean- 
field-type averaging procedure. The example shows that 
there may exist another mechanism: Nonlinearity via 
an entanglement pi| . Indeed a subsystem may start to 
evolve nonlinearly if its reduced density matrix evolves 
from a pure state to a mixture, and this happens when- 
ever the subsystem gets entangled with another one. 



S(p) = -(Tr(p)Tr(p 3 )) 



1/2 



(4) 



, is the "entropy" S3 given by Eq. (75) in (T|). The 
Lie-Nambu equation (|l|) is now equivalent to the matrix 
equation 



ip = 



Tr(p) 



1/2 r 



[h,p 2 



_Tr(p3) 

Its solution normalized by Tr p = 1 is 



1 

PW = 2 1 + ex P 



— iht 



3p 2 J 



p ■ cr exp 



iht 



3p 2 J 



(5) 



(6) 



For h — Eo~ 1 and p ■ cr = |cr 3 we find Ai = (1 + e)/2, 
A 2 = (1 - e)/2, and 



,1 _ ( cosu>(E,e)t 
A ~ \ —ismw(E,e)t 



-i siriuj(E, e)t 
cosui(E, e)t 



where ui(E,e) — 2E/y / l + 3s 2 . Notice that the vectors 
(j> A , k — 1,2, depend on Afc. In the linear case we can 
solve equations for orthogonal pure states and then form 
their convex combinations with coefficients A^ , which in 
no way affects the form of the pure states that form the 
mixture. In the nonlinear case the "pure state" compo- 
nents of the mixture do depend on the coefficients A& 
p9| . The dynamics of p(t) is nonlinear even though the 
Hamiltonian is given by the linear operator h. A possibil- 
ity of introducing nonlinearities without modifications of 
an algebra of observables is one of the important differ- 
ences between the Lie-Nambu formalism and nonlinear 
Schrodinger equations. For p 2 — 1/4 the density matrix 
is a projector and its dynamics is linear. 

Now consider a solution of ([!]) which at t — is a 
convex combination 



p(0) = pipi(0) + p 2 P2(0) 



(7) 



IV. POSITIVITY VS. COMPLETE POSITIVITY 

We have shown that self-adjoint Hilbert-Schmidt solu- 
tions of nonlinear Lie-Nambu equations are positive for 
1 7^ if they are positive at t — 0. Typically it is assumed 
that density matrices should be described by completely 
positive maps. A physical motivation behind complete 
positivity is the problem of extension of dynamics from 
subsystems to composite systems ]22] ] : One requires that 
a dynamics of a system described by p should allow to 
treat this system as a trivially embedded subsystem of 
a bigger one. The evolution of the bigger system should 
preserve positivity of its density matrix. If this is the 
case, and if for any t the density matrix of the composite 
system is positive independently of the dimension of the 
system we have added, then the dynamics of the original 
system is said to be completely positive. One additional 
technical assumption one makes is finite dimensionality 
of the system one adds. 

In the standard analysis of completely positive maps 
one assumes the maps are linear. But the dynamics we 
have discussed in this Letter is nonlinear. It can be shown 
that the definition of a nonlinear completely positive map 
introduced in mathematical literature p3| , p4| and applied 
to Hartree-type equations in |^6|,^5| is physically incor- 
rect. This problem is analyzed elsewhere p^Efl . From 
the point of view of this Letter it is sufficient to note that 
there exists a large class of nonliner triple-bracket equa- 
tions that do not lead to any problem with extension of 
dynamics from subsystems to composite systems. They 
satisfy all physical requirements typically associated with 
the notion of a completely positive map. Quite surpris- 
ingly, they are not in the form one takes as a departure 
point for the discussion of complete positivity of nonlin- 
ear maps in the mathematical literature. 

We gratefully acknowledge discussions we had on the 
subject with M. Kuna, P. Horodecki, G. A. Goldin and 
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K. Jones. The work of M. C. is a part of the Polish- 
Flemish project 007. 

V. APPENDICES 

A. Structure constants 

Consider a Hilbert space of vectors -0a where the ab- 
stract index a can be discrete, continuous, or composite 
p8[ . Denote the scalar product by (4>, tp) — uj aa ' 4> a -tp a ' ■ 
The tensor uj aa is in general a distribution whose inverse 
is I aa ' ■ By the inverse it is meant that 

w% = V"' (9) 
Lo aa 'l 0a , = 5 p a (10) 

where the <5's mean the Kronecker or Dirac deltas, or their 
products. In a Hamiltonian formulation of Schrodinger- 
type equations uj aa and I aa > play the roles of a symplec- 
tic form and a Poisson tensor, respectively The 
structure constants are 

Q a bc = 5p/ a <5 7 Q //3 7 ' - 5y a 5p a I-,i3> (11) 

Qabc = Ia0' I/3y' I-ya' ~ Iaj'Ipa'If/3' (12) 
n abc = _ [J «/3' W /37' W 7"' + ^j' ^a' ^0 _ ^3) 

Notice the spinor-type convention we use. Different equa- 
tions (nonrelativistic Schrodinger, positive-metric Dirac, 
off-shell Dirac, Bargmann-Wigner) correspond to differ- 
ent Hilbert spaces, w's and Ps but the form of the struc- 
ture constants is always the same. The indices are raised 
and lowered by metric tensors defined below. 

B. Metric tensors and Casimir invariants 

We define higher-order metric tensors by 

g ai...a n = 0J a i<cj a2a iUJ a3a 2 . . . u} a "-i<-2 UJ a n<-i (14) 
Gai ...a n I ai a' n 2 a ^ 3 a 2 ■ ■ ■ ^a n -i a ' n _o ^ a n a ' n _ 1 (*^) 

= 5o„...oi< (16) 

For n = lwe get just the symplectic form and the Pois- 
son tensor. For n — 2 we obtain the metric tensor on 
the Lie algebra — it is this tensor that lowers and raises 
the indices in the structure constants and (jl4|), ([l5|). For 
higher n's the tensors define higher order Casimir invari- 
ants of the Lie-Nambu bracket 

Cn(p) = g^-^Pa, ...Pa n = g ai ...a nP ai ■ ■ ■ (17) 

where p a = p aa ,. Let p a = ^a^a' be the 

Hilbert-Schmidt decomposition of p. The vectors 4>a are 
orthonormal (tJ aa (f>a4> l a , = rj kl — r] kk S kl ; rj kk = ±1 if 
the metric is not positive definite) which implies that 



C n {p) = E^iC^ )" A fe = Tr (P™)- The metric tensors 
are therefore an abstract index counterpart of trace. 
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